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Abstract 

A polarization analysis of the fine-structure intervals for the n = 17 Rydberg states of Mg and 
the n = 29 states of Si2+ is performed. The coefficients of all terms in the polarization expansion up 
to were computed using a semi-empirical single electron analysis combined with the relativistic 
all-order single-double method (MBPT-SD) which includes all single-double excitations from the 
Dirac-Fock wave functions to all orders of perturbation theory. The revised analysis yields dipole 
polarizabilities of ai = 35.04(3) a.u. for Mg+ and qi = 7.433(25) a.u. for Si'^"*', values only 
marginally larger than those obtained in a previous analysis (E. L. Snow and S. R. Lundeen (2007) 
Phys. Rev. A 75 062512, ihid (2008) 77 052501). The polarizabihties are used to make estimates 
of the multiplet strength for the resonant transition for both ions. The revised analysis did see 
significant changes in the slopes of the polarization plots. The dipole polarizabilities from the 
MBPT-SD calculation, namely 35.05(12) a.u. and 7.419(16) a.u., are within 0.3% of the revised 
experimental values. 

PACS numbers: 31.15.ap, 31.15.ag, 31.15.am, 31.15.V- 
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I. INTRODUCTION 



Resonant excitation Stark ionization spectroscopy (RESIS) [])] is a versatile and powerful 
method for studying Rydberg states of atoms and ions. One of the primary applications 
is the determination of deviations from the pure hydrogenic values of the binding energies. 
Polarization interactions between the core and the Rydberg electrons lead to the effective 
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This functional form has been applied to the analysis of the fine-structure spectrum of the 
Rydberg states of neutral Mg and Si^"*" resulting in precise estimates of the dipole polariz- 



abilities of the sodium-like Mg"*" and Si'^"'" ground states 



35.00(5) a.u. [6] and the 81^+ polarizability was 7.426(12) a 



u. 0. 



The Mg"*" polarizability was 



Analysis of the spectrum 



has also given information about the quadrupole polarizabilities. 

One area of uncertainty in the analysis is the contribution of the higher-order terms 
in the polarization expansion. Using theoretical estimates of C^ and CgL to constrain the 
analysis has proved essential in obtaining values of the quadrupole polarizability that are 
even remotely close with theoretical estimates j^, . However, some of the high order terms 
that contribute to Eq. ([T]) were omitted from the analysis of the experimental data. 

This limitation is rectified in the present work which uses two different theoretical tech- 
niques to determine values of all the terms in Eq. ([1]). One technique supplements the 
Hartree-Fock core potential with a semi-empirical polarization potential and effectively 
solves a one-electron Schrodinger equation to determine the excitation spectrum for the 
valence electron . The other method used is the relativistic all-order single-double 

method where all single and double excitations of the Dirac-Fock (DF) wave function are 
included to all orders of many-body perturbation theory (MBPT) |lOl . I 111 . Il2[ |. We note 
in passing that there has been a recent configuration interaction (CI) calculation of the 
polarizabilities of the Mg"*" and Si'^"*' ground states 1^. 

The current work has implications that go beyond the analysis of the RESIS experiments 
of the Lundeen group. One of the most active area in physics at present is the development 
of new atomic clocks jjased on groups of neutral atoms in optical lattices [l4, |l5, [iG] or 



single atomic ions 



0. 



These clocks have the potential to exceed the precision of 
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the existing cesium microwave standard 



For many of these clocks the single largest 



19 



20 



21 



22 



23| 



source of systematic error is the black-body radiation shift (BBR) 
The BBR shift to first order is proportional to the difference in polarizabilities of the two 
states involved in the clock transition. Many estimates of the relevant polarizabilities are 



determined by theoretical calculations 2J, |25|, |26|] . Comparisons of existing techniques to 



calculate polarizabilities with high quality experiments will ultimately help constrain the 
uncertainties associated with the BBR shift. 

II. THE POLARIZATION EXPANSION 



In this section the definitions of the various terms in the polarizationpotential are given 
following the analysis of Drachman . The notation of Lundeen [l|, |5j is adopted. 
The leading term, C4 is half the size of the static dipole polarizability, 



The dipole polarizability is defined as 
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where fgn is the absorption oscillator strength for a dipole transition from state g to state 
n. The absorption oscillator strength for a multi-pole transition from g ^ n, with an energy 



difference of AE^g = Eg — En, is defined as 
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In this expression, Lg is the orbital angular momentum of the initial state while k is the 
polarity of the transition. In a J-representation, the oscillator strength becomes 
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The next term, Cg, is composed of two separate terms 
The quadrupole polarizability, 02 is computed as 
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The second term in Eq. ([6]) is the non-adiabatic dipole polarizabihty. It is defined as 



The r term, C7 also comes in two parts, namely 
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The 7i is a higher-order non-adiabatic term 
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while q is the charge on the core. The dipole-dipole-quadrupole polarizabihty, ai p ar ises 
from third order in perturbation theory. It is derived from the matrix element |2|. Isl. |27| 

aii2 {^g;0\V''^\tpna]La) 
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where = C'^(f)-C^(R)r'^/i?'^"'"^. The sum of the multipole orders must obey ki + k2 + k3 = 
4. Quite a few terms contribute to Cg 

as- P2- aiPi + ami + 
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The octupole polarizabihty, as is computed as 

«3 = E 
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The comes from the non-adiabatic part of the quadrupole polarizabihty, it is 
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The fourth-order term, aim is related to the hyper-polarizability [28|, l29|]. It is defined as 
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The final term, C^l is non-adiabatic in origin and defined 
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III. STRUCTURE MODELS FOR Mg+ AND Si^+ 



A. Semi-empirical method 



The semi-empirical wave functions and transition operator expectation values were com- 

y, 
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in a large mixed 



puted by diagonalizing the semi-empirical Hamiltonian 
Laguerre type orbital (LTO) and Slater type orbital (STO) basis set [30|]. We first discuss 
and then mention Mg"*". 
The initial step was to perform a Hartree-Fock (HF) calculation to define the core. The 
present calculation can be regarded as HF plus core polarization (HFCP). The calculation 
of the Si^^ ground state was done in a STO basis 3J]. The core wave functions were then 
frozen, giving the working Hamiltonian for the valence electron 

1. 



H = + V.Jr) + K.c(r) + VJr) 
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The direct and exchange interactions, Vjir and Vexc, of the valence electron with the HF core 
were calculated exactly. The ^-dependent polarization potential, Vp, was semi-empirical in 
nature with the functional form 
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The coefficient, ad is the static dipole polarizability of the core and gj{r) = 1 — exp(— r^/p^) 
is a cutoff function designed to make the polarization potential finite at the origin. The cutoff 
parameters, pi were tuned to reproduce the binding energies of the ns ground state and the 
np, nd and nf excited states. The dipole polarizability for Si^+ was chosen as = 0.1624 
35]. The cutoff parameters for £ = ^ 3 were 0.7473, 0.8200, 1.022 and 0.900 qq 



a.u. 



30 



respectively. The parameters for £ > 3 were set to ps. The energies of the states with ^ > 1 
were tuned to the statistical average of their respective spin-orbit doublets. The Hamiltonian 
was diagonalized in a very large orbital basis with about 50 Laguerre type orbitals for each 
£-value. The oscillator strengths (and other multi-pole expectation values) were computed 



with operators that included polarization corrections [30|, iMl, l36|, l37|, l38j. The quadrupole 
core polarizability was chosen as 0.1021 a.u. [3^] while the octupole polarizability was set 
to zero. The cutoff parameter for the polarization correction to the transition operator was 
fixed at 0.864 ao (the average of po. Pi, p2 and ps). 
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It is worth emphasizing that model potential is based on a realistic wave function and the 
direct and exchange interactions with the core were computed without approximation from 
the HF wave function. Only the core polarization potential is described with an empirical 
potential. 

The overall methodology of the Mg"*" calculation is the same as that for Si^+ and many of 
the details have been given previously j^]. The core dipole polarizabilities were ad = 0.4814 
a.u. 0, and Oq = 0.5183 a.u. for Mg^"*" 30, 35|. The octupole polarizability was set 
to zero. The Mg^+ cutoff parameters for £ = ^ 3 were 1.1795, 1.302, 1.442, and 1.520 
Oo respectively. The cutoff parameter for evaluation of transition multipole matrix elements 
was 1.361 flo- 

The HFCP calculations of the polarizabilities utilized the list of multipole matrix elements 
and energies resulting from the diagonalization of the effective Hamiltonian. These were 
directly used in the evaluation of the polarizability sum rules. 



B. The all-order method 

In the relativistic all-order method including single, double, and valence triple excitations, 
the wave function is represented as an expansion 
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where is the lowest-order atomic state wave function, which is taken to be the frozen-core 
DP wave function of a state v in our calculations. In second quantization, the lowest-order 
atomic state function is written as 

1$^) = 4|0c), 

where \0c) represents the DF wave function of the closed core. In Eq. ([inD, 4 and are 
creation and annihilation operators, respectively. Indices at the beginning of the alphabet, 
a, b, ■ ■ ■ , refer to occupied core states, those in the middle of the alphabet m, n, ■ ■ ■ , refer to 
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excited states, and index v designates the valence orbital. We refer to pma, Pmv as single core 
and valence excitation coefficients and to pmnab and pmnva as double core and valence excita- 
tion coefficients, respectively. The quantities Pmnwab are valence triple excitation coefficients 



and are included perturbatively where necessary as described in Ref. 

To derive the equations for the excitation coefficients, the wave function given by 
Eq. (fT9|) . is substituted into the many-body Schrodinger equation 

H\^,) = E\^,), (20) 



391] . This can be expressed 



where the Hamiltonian H is the relativistic no-pair Hamiltonian 
in second quantization as 

= ^ ei : alai : +^^gijki ■ alajazflfc :, (21) 

i ijkl 

where is the DF energy for the state i, gijki are the two-body Coulomb integrals, and 
: : indicates normal order of the operators with respect to the closed core. In the no-pair 
Hamiltonian, the contributions from negative-energy (positron) states are omitted. 

The resulting all-order equations for the excitation coefficients pma, Pmv, Pmnab, and pmnva 
are solved iteratively with a finite basis set, and the correlation energy is used as a conver- 
gence parameter. As a result, the series of correlation correction terms included in the SD 
(or SDpT) approach are included to all orders of many-body perturbation theory (MBPT) 
as an additional MBPT order is picked up at each iteration. The basis set is defined in a 
spherical cavity on a non-linear grid and consists of single-particle basis states which are lin- 



ear combinations of 5-splines 40|] . The contribution from the Breit interaction is negligible 
for all matrix elements considered in this work. 

The matrix element of any one-body operator Z in the all-order method is obtained as 

The numerator of the resulting expression consists of the sum of the DF matrix element 
and twenty other terms Z^^\ k = a ■ ■ - t. These terms are linear or quadratic functions 
of the excitation coefficients pma, Pmv, Pmnab, and Pmnva- More details on the SD and 
SDpT methods and their applications can be found in Refs. JjJ, |41[ . We find that the 
contribution of triple excitations is small for the atomic properties considered in this work. 
So the SD approximation is used for most transitions. 



The S-spline basis used in the calculations included = 50 basis orbitals for each angular 
momentum within a cavity radius of Rq = 100 ao for Mg"*" and Rq = 80 for Si^"^. Such 
large cavities are needed to fit highly-excited states such as 8h needed for the 3d octupole 
polarizability calculations. The single-double (SD) all-order method yielded results for the 

primary ns — npj electric-dipole matrix elements of alkali- metal atoms that are in agreement 

i 

with experiment to 0.1%-0.5% [11]. We refer to the results obtained with this method as 
MBPT-SD in the subsequent text and Tables. 

Since the all-order calculations are carried out with a finite basis set, the sums given by 
Eqs. ([3]) - ( |T3l) run up to the number of the basis set orbitals (A^ = 50) for each partial wave. 
For consistency, the same 5-spline basis is used in all calculations of the same system (e.g. 
Mg+ or Si3+). 

The calculation of the polarizabilities for the MBPT-SD uses slightly different pro- 
cedures to include different parts of the polarizability sum rules. The all-order matrix 
elements were combined with the experimental energies for excited states with n < 6 
for jS = ns,npi/2,np3/2,nd3/2,nd5/2, n < 7 for /3 = nf5/2, nf j/2, and n < 8 for /3 = 
ngi/2,ngg/2,nhQ/2,nhii/2- The remaining matrix elements and energies were calculated in 
the DF approximation, with the exception of the 3s dipole polarizability, where the remain- 
ing matrix elements were calculated using random-phase approximation (RPA) [42] for the 
purpose of error evaluation. These remainder contributions are small for dipole polariz- 
abilities (0.2-5%) but increase in relative size for the quadrupole (0.3-10%) and octupole 
(4-20%) polarizabilities. An extra correction was introduced to the remainder contribution 
for octupole polarizabilities. First, the accuracy of the DF calculations was estimated from a 
comparison of the DF and all-order results for the few first terms. Then, these estimates were 
used to adjust the remainder. The improvement of the DF results for states with higher n 
was also taken into account. The size of this extra correction ranged from 0.9% to 6% of the 
tail contributions as the accuracy of the DF approximation for these highly-excited states is 
rather high. The net effect of this scaling was usually to reduce the octupole polarizabilities 
by an amount of about 0.5-1.5%. 

The core contribution was calculated in the RPA [35,] with the exception of the dipole 
polarizability for the Mg^+ core. In this case the polarizability of = 0.4814 a.u. was 



taken from a pseudo-natural orbital CI type calculation 
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30|. A small correction for 



the dipole polarizability that compensates for excitations from the core to occupied valence 
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states was also determined using RPA matrix elements and DF energies. The relative impact 
of the core polarizability was at least a factor of two smaller for the quadrupole polarizability. 



IV. GROUND AND EXCITED PROPERTIES 



A. The energy levels 

The binding energies of the low-lying states of the Mg"^ and Si^'^ are tabulated and 
compared with experiment in Table [H The agreement between the HFCP energies and the 
experimental energies is generally of order 10~^ Hartree. When the pe cutoff parameters are 
tuned to the lowest state of each symmetry the tendency is for higher states of the same 
symmetry to be slightly under-bound. The MBPT-SD binding energies generally agree with 
experiment to better than 10~^ Hartree. The MBPT-SD binding energies do not suffer any 
systematic tendency to either underbind or overbind as n increases. 



B. Line strengths 



Table HDlists the line strengths for the resonant transitions of Na, Mg"*", AP+ and Si^"*". All 
line strengths here and in the text below are given in a.u.. The HFCP values for sodium are 
from calculations previously reported in Ref. 58|] while the values for AP"*" were taken from 



a calculation very similar in style and execution to the present calculations [59|. The MBPT- 
SD line strengths for Na and AP+ were taken from Ref. 10|. Values from the extensive 
tabulation of dipole line strengths using a 5-spline non-orthogonal configuration interaction 
with the Breit interaction (BSR-CI) [60i] are also listed. The HFCP line strengths were 



computed from a common multiplet strength by multiplying by the appropriate recoupling 



coefficients 



61|. 



The comparisons for the resonant 3s —* 3p transition reveal that the HFCP line strengths 
are the smallest, the BSR-CI line strengths are the largest and the MBPT-SD line strengths 
are intermediate between these two calculations. The MBPT-SD dipole strengths are closer 
to HFCP for Na, Mg+ and A\'^+ and about half-way between HFCP and BSR-CI for 
The total variation between the three different calculations is about 1%. The most precise 



experiments performed on 



on sodium itself 
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le Na-like iso-electronic series of atoms(ions) are those performed 



481]. The experimental line strengths for sodium are in better 



agreement with the MBPT-SD and HFCP hne strengths than they are with the BSR-CI 
hne strengths. 

There have been two precision measurements of the 3s —* 3p transition rate for Mg"*". 



The experiment of Ansbacher et al. 



49j gave shghtly larger hne strengths which agree best 



with the BSR-CI values. However, the most recent trapped ion experiment [50;] gave a 
3si/2 — ^ 3p3/2 line strength of 11.24(6) that is in better agreement with the HFCP/MBPT- 
SD line strengths. 

Table IIIII lists the line strengths for a number of other dipole transitions for Mg+ and 
Si^+. The line strengths for the quadrupole 3s nd transitions are also listed due to their 
importance in the determination of the quadrupole polarizabilities. 

The 3p 3d transition is the strongest transition emanating from the 3p level. The 
comparison between the three calculations exhibits a pattern similar to that of the resonant 
transition. The HFCP line strengths are smallest, the BSR-CI line strengths are the largest, 
and the MBPT-SD line strengths lie somewhere between these two calculations. 

The astrophysically important Mg~^ 3s — > Ap transition has a very small dipole strength. 
It is close to the Cooper minimum 62|] in the 3s — np matrix elements and therefore is more 
sensitive to the slightly different energies between the spin-orbit doublet. This caused the 
ratio of line strengths for the 4pi/2 and 4p3/2 transitions to deviate from the expected value 
of 2. The MBPT-SD branc hing ratio of 1.76 agrees with the recent experimental values of 
1.74(6) y and 1.82(8) The HFCP muhiplet strength of 0.00752 and the MBPT- 

SD multiplet strength of 0.00721 are about 5-10% smaller than the recent experimental 
estimates of 0.00793(26) [esj and 0.00775(50) ^M\. 

There is also a deviation from the ratio of 2 for the 3s Api/2.^/2 transitions of Si'^"^. 
However, in this case the deviation is smaller. Ratios of line strengths for the stronger 
transitions are much closer to values expected from purely angular recoupling considerations. 
The 3p3/2 : 3j9i/2 ratio for Si'^^ was 2.002. The 3p 4s transition ratio has a slight deviation 
from 2, the MBPT-SD calculations giving 2.015 for Mg+ and 2.006 for Si^+ (the BSR-CI 
ratios are similar). 

The better than 0.5% agreement between the model potential and MBPT-SD line 
strengths for strong transitions is consistent with previous comparisons. The general level 
of agreement between calculations with a semi-empirical core potential and more sophis- 
ticated ab-initio approaches for properties such as oscillator strengths, polarizabilities and 
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dispersion coefficients has generally been very good [30|, l65|, l66|, |67|]. There was a tendency 
for the agreement between the HFCP and MBPT-SD line strengths to degrade slightly from 
Mg"*" and Si^"*". This is probably due to the increased importance of relativistic effects as the 
nuclear charge increases. 



C. Polarizabilities 

The polarizabilities of the 3s, 3p and 3d levels of Mg^ and Si^"*" are listed in Table llVi 
Tensor polarizabilities are also determined for the 3p and 3d levels. Definitions of the tensor 



68| and 



)olarizability, ai,2jj, in terms of oscillator strength sum rules can be found in Refs. 

Table |V] gives a short breakdown of the contributions of different terms to the dipole 
polarizability while Table |Vl] gives the breakdown for the quadrupole polarizability. The 
3s ep{d) contribution represents anything over n = 6 and can be regarded as a mix of 
some higher discrete states as well as the pseudo-continuum. Polarizabilities for the Mg"*" 
and Si^"*" ground states from other sources are also listed in Table |V] and IVII The HFCP 



Mg"*" polarizability is marginally smaller than that reported previously 25|] since the present 
evaluation includes a small core- valence correction. 

The very good agreement between the HFCP and MBPT-SD polarizabilities is a notable 
feature of Table HVl None of the static polarizabilities differ by more than 0.5% with the 
exception being the a2 of the Mg"*" 3d state. Here the difference is caused by the very 
small Ai?3(i_4s energy difference which is sensitive to small errors in the HFCP energies. 
The relative difference between some of the tensor polarizabilities is larger, but this is due 
to cancellations between the component sum rules that are combined to give the tensor 
polarizability. 

A recent CI calculation of the Mg+ and Si^+ ground state dipole polarizabilities [l 
gave polarizabilities that were 1-2% larger than the HFCP/MBPT-SD polarizabilities. The 
more recent relativistic coupled-cluster (RCC) calculation [52,] gave polarizabilities that were 
compatible with the present values. 

The dipole polarizabilities for both Mg+ and Sr^^ are dominated by the resonant oscillator 
strength. For Mg+ one finds that 98.3% of ai arises from the 3s —>■ 3p transition. For 
Si^+ the contribution is smaller but still substantial at 96.7%. The non-adiabatic dipole 
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polarizabilities are even more dominated by the contribution from the resonant transition. 
One finds that 99.9% of jSi and 99.99% of 71 for Mg"*" come from this transition. The 
proportions for the Si^"*" /?i and 71 are 99.6% and 99.92% respectively. 

The quadrupole polarizabihties are also dominated by a single transition. Table IVTl shows 
that the 3s 3d excitation constitutes at least 95% of 02 for both Mg"*" and Si^+. 

The calculation of the aiu and aim polarizabilities was a composite calculation using 
both MBPT-SD and HFCP matrix elements. The HFCP calculation automatically gener- 
ates a file containing matrix elements between every state included in the basis. The more 
computationally intensive MBPT-SD calculation was used for the largest and most impor- 
tant matrix elements. The HFCP matrix elements for the 3s — 3p, 3p — > 3d, 3s 3d 
and 3p 4s transitions were replaced by J weighted averages of the equivalent MBPT-SD 
matrix elements. This procedure combines the higher accuracy of the MBPT-SD calculation 
with the computational convenience of the HFCP calculation. The justification for this pro- 
cedure is that the predominant contribution to the polarizability comes from the low-lying 
transitions. The resulting polarizabilities are listed in Table IVlIl The biggest change in an2 
and aim resulting from using the composite matrix element list was less than 0.3%. The 
aii2 and aim polarizabilities did not allow for contributions from the core. The impact of 
the core will be small due to large energy difference involving core excitations. The relative 
effect of the core for aii2 and aim can be expected to be about as large as the core effect in 
the ground state ai and a2 since there are core excitations that contribute to with only one 
core energy in the energy denominator. For example, consider the aii2 excitation sequence 
of 2/3s ^S^ 2p^3s3d ^F" 2p^3d ^D^ ^ 2/3s ^S^ The numerical procedures used to 
generate the aii2 and aim polarizabilities were validated for He^. A calculation of the He^ 
excitation spectrum was performed and the resulting lists of reduced matrix elements were 
entered into the polarizability programs. All the coefficients given by Drachman TOj were 
reproduced. 

The polarizabilities in Table IVlTl from the composite matrix element list could be regarded 
as the recommended set of polarizabilities. The MBPT-SD matrix elements are used for the 
dominant low-lying transitions. The HFCP matrix elements are more accurate than the 
RPA/DF matrix elements used for the 3s ep{d) remainders. The only difference between 
the Table IVIII and MBPT-SD polarizabilities occurs for a^. A relatively large part of of 
as comes from the higher excited states and the continuum. Accumulating a lot of small 
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contributions is tedious for the computationally expensive MBPT-SD, so this is done with 
the less accurate DF approach. In this case the HFCP polarizability is to be preferred. It 
should be noted that the octupole polarizability is of minor importance in the subsequent 
analysis. 



D. Error assessment 



Making an a-priori assessment of the accuracy of the HFCP polarizabilities is problematic 
since they are semi-empirical in nature. The error assessment for the MBPT-SD proceeds by 
assuming that the total contribution of fourth- and higher-order terms omitted by the SD all- 
order method does not exceed the contribution of already included fourth- and higher-order 
terms. Thus, the uncertainty of the SD matrix elements is estimated to be the difference 
between the SD all-order calculations and third-order results. 

This procedure was applied to the S^s-spi/i ^^^^ strength of sodium yielding an uncertainty 
SS3s-3pi^2 ~ 0.092. This uncertainty exceeds the difference between the SD line strength 
of 12.47 |lO| and recent high precision experiments which give 12.412(16) 4J, |45|, and 
12.435(41) |46|]. A similar situation applies for the S's^-spg^j line strength. 

A detailed first principles evaluation of the uncertainty of the Si^^ static dipole polar- 
izability has been done and the uncertainty budget is itemized in Table El In this case, 
the difference between the SD line strength and third order line strength for the resonance 
transition was 0.082%. The uncertainties in the remaining (n = 4 — 6) discrete transitions 
were of similar size. Uncertainties in the energies used in the oscillator strength sum rule 
can be regarded as insignificant since experimental energies were used. To estimate the 
accuracy of the remainder of the valence sum, the (n = 4 — 6) calculation was repeated 
using RPA matrix elements and DF energies. The difference of 3% between the MBPT-SD 
and DF/RPA values was assessed to be the uncertainty in the ep remainder. The good 
agreement between the HFCP and DF/RPA for the non- resonant valence contribution gives 
additional evidence that the uncertainty estimate is realistic. 

The core dipole polarizability calculated in the RPA is known to underestimate the actual 
core polarizability. For neon, the RPA gives ai = 2.38 a.u. 3^ which is 11% smaller than 



the experimental value of 2.669 a.u. 7l|]. For Na"^, the RPA gives 0.9457 a.u. 35] while 
experiment gives 1.0015(15) a.u. 72|. The pseudo-natural orbital approach used for Mg^+ 
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Q and ai = 0.9947 a.u. for Na+ Q- 



j;ave ai = 2.67 a.u. for Ne [73|] and ai = 0.9947 a.u. for Na+ [Tj. The uncertainty in the 



quadrupo 



neon 
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e core polarizabihty is based on comparisons with coupled cluster calculations for 



75| . The RPA value of 6.423 a.u. is about 12% smaller than the coupled cluster 



values of 7.525 a.u. [75] and 7.525 a.u. [74]. The relative uncertainties are 5ai(Mg ) = 2%, 



5ai(Si'^^) = 5%, 6a2(Mg^^) = 12%, and 6a2{Si'^'^) = 12%. The core-valence correction was 
assigned an uncertainty of 20% based on differences between DF and RPA matrix elements. 
The RPA error estimates are likely to be very conservative since the uncertainty in the RPA 
polarizabilities is expected to decrease as the nuclear charge increases. 

Combining the uncertainties in the valence and core polarizabilities for Si^~^ gives a final 
uncertainty of 0.16 a.u. (or 0.22%) in the MBPT-SD ai. 

The uncertainty in the Si^"*" 0^2 listed in Table I VI I was evaluated with a process that was 
similar to the dipole polarizabihty. The difference between the SD line strength and third 
order line strength for the 3s —> 3^5/2 transition was 0.064% (the relative uncertainty was 
almost the same for the transition to the 3^3/2 state). This uncertainty is slightly smaller 
than that for the resonant dipole transition. This was expected since the 3d electron is 
further away from the nucleus than the 3p electron and therefore correlation-polarization 
corrections have less importance. Rather than do a computationally expensive analysis, the 
relative uncertainties in the {nd + ed) remainders were conservatively assigned to be same 
as for the dipole transitions. The final uncertainty was 6a2 = 0.03 a.u.. 

The relative uncertainties in the Mg"*" polarizabilities are set in the same way as Si^"*". 
The difference between the third-order and all-order dipole line strengths for the resonance 
transition was 0.3%. The relative differences were larger for the n = 4 — 6 transitions due to 
their small size. For example, the third-order/all-order comparison for the S^g^p multiplet 
strength gave 5%. This is consistent with the difference with the experimental multiplet 
multiplet strength. The uncertainties were slightly smaller for the slightly larger 5p and 6p 
transitions. However, the net contribution to the uncertainty was miniscule since the line 
strengths were so small. The 3s —>■ ep uncertainty of 5% was based on differences between 
the HFCP and DF/RPA matrix elements. 

The uncertainties in the Mg"^ 02 polarizabihty are listed in Table I VI I and n = 3 — 6 
transitions were derived from the third-order /all-order comparison. The relative uncertainty 
in the 3s — > 3d transition was 0.22%. The very good agreement between the HFCP and 
MBPT-SD values for these terms is further supportive of a small uncertainty for the n = 3—6 
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transitions. The 7% uncertainty in the 3d ed remainder was based on the differences 
between the MBPT-SD and DF matrix elements. 

The relative uncertainties in the octupole polarizabilities listed in Table IVIII were set to 
the uncertainties in the quadrupole polarizabilities. The nf orbitals are further away from 
the core than the 3c/ orbitals and so the 02 uncertainty serves as a convenient overestimate. 

The uncertainties in the higher-order polarizabilities /?2 and 71 listed in Table IVIII 
were taken to be the uncertainties in the resonant line strengths. The higher powers in the 
energy denominator means other transition make a negligible contribution. 

The uncertainties in aii2 and aim were derived from the uncertainties in the reduced 
matrix elements. The relative uncertainties for the most important 3s 3]9, 3]9 — > 3d and 
3s 3d matrix elements were simply added to give relative uncertainties for valence part 
of aii2 and aim. The relative uncertainty resulting from the omission of core excitations 
was taken as the ratio of the core to total dipole polarizability and added to the aii2 and 
aim uncertainties. 

The uncertainties in Ce, C7, Cg and C^l were determined by combining the uncertainties 
of the constituent polarizabilities. The most important of these parameters is the expected 
slope of the polarization plot, i.e. 5Cq = 5a2/2 + 35/5i. For Si^^ we get 6Ce = 0.015 + 0.027 = 
0.042. For Mg+ the uncertainty was dC^ = 1.2. 

V. POLARIZATION ANALYSIS OF RYDBERG STATES 
A. The polarization interaction 

The various polarizabilities needed for the polarization analysis are listed in Table IVIII 
The Cj, Cs and Csl values were used to make corrections to the experimental energy in- 
tervals. The C4 value was used in computing the second-order energy shift. The transition 
matrix elements used in this calculation represent a synthesis of the HFCP and MBPT-SD 
calculations. 

Estimates of C7 and Csl were previously made by Snow and Lundeen using MBPT-SD 
transition amplitudes for the lowest lying transitions. These earlier estimates are within a 
few percent of the present more sophisticated analysis. The Snow and Lundeen values for 
C7 were —1684(9) a.u. for Mg"*" and —122(9) a.u. for Si^"*". They are a few percent smaller 
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than those hsted in Table IVIII due to the omission of higher excitations from the sum rule. 
The Snow and Lundeen values for C^l were 1170(12) a.u. for Mg+ and 60.5 a.u. for Si'^+. 

One aspect of Table IVIII that is relevant to the interpretation of experiments is the 
importance of the non-adiabatic dipole polarizabilities. Consider Mg"*" for example. The 
respective contributions to Cg are 78.05 a.u. from a2 and —318.0 a.u. from — 6/?2- Similarly, 
one finds that the 71 term of —1.6 x 324.7 makes up 30% of the final Cj value of —1727 a.u. 
And finally, one finds that the Cg value of 10672 is largely due to the 3671 contribution of 
11689 a.u.. The degree of importance of the non-adiabatic terms scarcely diminishes for the 
Si^^ ion. 

Table IVIIII gives the energy shifts to the n = 17 levels of Mg"^ and the n = 29 levels of 
Si'^+ using the values in Table IVIII The energy shifts need (r~") expectation values which 
were evaluated using the formulae of Bockasten 



B. The polarization plot 

Polarizabilities can be extracted from experimental data by using a polarization plot. 
This is based on a similar procedure that is used to determine the ionization limits of atoms 



771]. The notations -B4 and -Be, (instead of C4 and Ce) are used to represent the polarization 
parameters extracted from the polarization plot. This is to clearly distinguish them from 
polarization parameters coming from atomic structure calculation. Assuming the dominant 
terms leading to departures from hydrogenic energies are the i?4 and terms, one can write 

^4 + i?6-rj-7f . (23) 



A(r-4) " "A(r-4) 
In this expression, Ai? is the energy difference between two states of the same n but different 
L, while A(r^^) and A(r^^) are simply the differences in the radial expectations of the two 
states. 

There are other corrections that can result in Eq. fl23|) departing from a purely linear 
form. These are relativistic energy shifts. Stark shifts due to a residual electric field, and 
polarization shifts due to the C-j, Cg (and possibly higher-order) terms of Eq. ([T]). The 
energy difference between the (n, L) and (n, L') states can be written 

AE = AE4 + ^Eq + AE7 + AE-g + AEgi 

+ AE,el + AEsec + AE,s , (24) 
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where arises from the polarization terms of order (r "■). 

Dividing through by A(r~^) and replacing Ai^g by B^A^r"^) gives 

-D4 + -D6- 



A(r-4) " "A(r-4) A(r-4) 

AErcl + AEsec + AEss 

+ A{^) • (2^) 

The influence of the Stark shifts, relativistic shifts, and second-order polarization correc- 
tion can be incorporated into the polarization plot by simply subtracting the energy shifts. 
The corrected energy shift, AE^i, is defined as 



AEci AEohs AEj-ci + AEscc + AE^s 



(26) 



A(r-4) A(r-4) A{r-^) 

An approximate expression is used for the relativistic energy correction. This is taken from 
the result 

E,,^-^('-±]. (27) 

2n^ \j + l/2 An) ^ ' 

The correction due to second-order effects, AE^pr, uses the results of Drake and Swainson 



The Stark shift corrections use the Stark shift rates from Snow and Lundeen 5, 6] and 
the deduced electric field. The energy corrections due to relativistic and polarization effects 
for the states of Mg+ and the Si"^^ for which RESIS data existed are listed in Table IVIIIi 

The second corrected energy is defined by further subtracting the polarization shifts, 
AE7, AEg and AEgL, 

AE,2 AE,i AEj + AEs + AEsl 



A(r-4) A(r-4) A{r-^) 

C. Mg+ 



(28) 



The energy splitting between adjacent L Rydberg levels is dominated by the C4 term. 
The next biggest term is the AEq term which is 3% of A£'4 for the (17,6)-(17,7) interval. 
The AEgL correction is larger than AE^. The relative impact of the higher-order corrections 
diminishes as L increases. 

The revised analysis of the RESIS energy intervals for Mg"*" was performed by subtracting 
the Ai^ci and Aii^c2 energy corrections itemized in Table IVIIII from the observed energy 
splittings. This represents a refinement over the previous analysis by Snow and Lundeen 
6| in a couple of respects. First, Snow and Lundeen did not include the Cg term since the 
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necessary polarizability information simply was not available. Their evaluation of aiu only 
included the 3p and 3d states in the intermediate sums. The truncation of the sums in the 
aii2 calculation was justified as the correction to aii2 from a more complete evaluation was 
only a few percent. The impact of the AEg shift is more substantial. Table IVlnl shows the 
relative size of AE^ with respect to AEj ranging from 35% to 11%. 

Figures [1] shows the polarization plot for Mg"*". Linear regression was applied to the four 
data points with A(r~^)/A(r~^) < 0.002. The (17,6)-(17,7) interval was omitted from the 
fit because the influence of AEj s^ql and Ai^sec amount to just over 50% of AEq. Visual 
examination of Figure [1] shows this data point lies a significant distance away from the line 
of best fit obtained from the four remaining points. The linear regression gave an intercept 
of i?4 = 17.522(7) a.u. and a slope of Bq = —251.2(79) a.u.. The quoted uncertainties are 
the statistical uncertainties from the linear regression fit. 

The new value of the dipole polarizability derived from the polarization plot intercept 
was 35.044 a.u.. This is marginally larger than the polarizability of 35.00(5) a.u. given in 
the original Snow and Lundeen analysis The present ai is larger because the additional 
corrections in the AEc2 energies lead to a steeper polarization plot. 

The slope of Bq = —251.2(79) is slightly steeper than the Table IVlTl recommended Cq 
of —240.1(12). Using the slope of —251.2 in conjunction with a (3i = 106.0 a.u. gives a 
quadrupole polarizability of 0:2 = —502.4 + 636.3 = 133.9 a.u.. This is about 90 a.u. smaller 
than the polarizability of 222(54) a.u.. given by Snow and Lundeen j^l. However it is only 
22 a.u. smaller than the theoretical polarizabilities of 156.1 a.u.. The uncertainty in the 
derived quadrupole polarizability would be (2 x 7.9 + 6 x 0.3) = 17.6. The RESIS and 
theoretical values are slightly outside their respective combined error estimates. However, 
the uncertainty estimate used for Cq is purely statistical in nature and does not allow higher 
order corrections to Eq. (!23|) . This point is discussed in more detail later. 

The relatively large change in 02 from 222(54) to 134(18) a.u. was caused by the inclusion 
of AEg. There is a near cancellation between some of the AEj and AEgi energy corrections. 
Hence the inclusion of the AEg energy correction has a relatively large impact. For example, 
the sum of AE-i and AE^l for the (17,7)-(17,8) interval was 1.473 MHz. The AEg correction 
was 0.851 MHz. 

The derived dipole polarizability and value of Bq are not sensitive to small changes 
in the Cn values used for the corrections. An analysis using alternate C„ values derived 
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from the uncertainties detailed in Table IVIII was performed. This resulted in an additional 
uncertainty of 0.0004 a.u. in and an additional uncertainty of 1.6 in Bq. These additional 
uncertainties were sufficiently small to ignore in subsequent analysis. 



17.6 

^ 17.4 
< 

W 17.2 

<l 

17.0 
16.8 

0.001 0.002 0.003 

A<r"^)/A<r"^) 

FIG. 1: The polarization plot of the fine-structure intervals of Mg for the n = 17 Rydberg levels. 
The AEci intervals are corrected for relativistic, second-order and Stark shifts. The AEc2 intervals 
account for and (r^®) shifts. The linear regression for the A£'c2 plot did not include the last 

point. 





0.002 0.004 0.006 0.008 

A<r"^>/A<r"^> 

FIG. 2: The polarization plot of the fine-structure intervals of Si^"*" for the n = 29 Rydberg levels. 
The AEci intervals are corrected for relativistic, second-order and Stark shifts. The AEc2 intervals 
account for (r^^) and (r~®) shifts. The linear regression for the AEc2 plot did not include the last 
two points. 
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D. Si3+ 



The polarization plot for Si^^ is shown in Figure [2l The most notable feature is the 
large difference between the Aii^ci and AEc2 data-sets. The other notable feature is the 
pronounced deviation from linear of the plot. 

Examination of Table IVIIII for the (29,8)-(29,9) interval shows that the net Ai?7 8,8L 
correction is very close in magnitude to the AEq energy correction. The AEy^s,sl correction 
is still more than 50% of the AEf, correction for the (29,10)-(29,11) interval. The polarization 



series is an asymptotic series 



m 



79| and is not absolutely convergent as n increases. As 



mentioned by Drachman [70], a condition for the usefulness of the polarization series is 
that the AEj^s^sl corrections should be significantly smaller than the AEq corrections. This 
condition is not satisfied for the first two intervals and leads to the noticeable curvature in 
the plot of the AEc2 data points. 

The resolution to this problem would be to increase the L values at which the intervals 
are measured. But Stark shift corrections become increasingly important at high L. The 
Stark shift corrections are significant for the (29,11)-(29,14) interval. 

A line of best fit was drawn using the four data points with A{r~^) / A{r~^) < 0.004. The 
linear regression gave an intercept of -B4 = 3.7163(32) and a slope of Bq = —30.96(134). 
The intercept translates to a polarizability of 7.433 a.u.. To put this in perspective, the 
polarizability originally deduced from the RESIS experiment was 7.408(11) ^. A later 
analysis which included the C7 and Csl potentials gave 7.426(12) a.u. j^. There has been 
a steady increase in the derived dipole polarizability as more higher-order terms in the 
polarization series are incorporated into the analysis. 

The polarization plot Bq of —30.96(134) was about 10% larger in magnitude than the 
MBPT-SD value of —27.06(5). This value of Bq results in a quadrupole polarizability of 
a2 = (-2 X 30.96 + 6 X 11.04) = 4.34 a.u. which is 60% smaller than the HFCP and MBPT- 
SD polarizabilities. The uncertainty of (2 x 1.34 + 6 x 0.006) = 3.0 a.u. is too small to allow 
consistency with the theoretical values. 
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E. An alternate perspective 



The analysis so far can be regarded as a standard polarization analysis but with additional 
refinements due to improved knowledge about the higher-order terms in the polarization 
series. However, it is worthwhile to examine the analysis from a different perspective. 

The comparison between first principles theory and the RESIS experiment has resulted 
in agreement to better than 1% for dipole polarizabilities. The quality of the agreement 
for the quadrupole polarizability is not nearly so good. But can the analysis of the RESIS 
experiment be expected to yield quadrupole polarizabilities that are a serious test of calcu- 
lation? The quadrupole polarizability is derived from the slope of the polarization plot. But 
the higher-order polarization corrections and Stark shifts result in energy corrections that 
amount to between 30-100% of the raw Cq energy shift. And it must be recalled that the 
polarization series itself is an asymptotic series 79] so there are uncertainties about the size 
of omitted terms. 

One way forward is to use the dipole polarizability comparison as a guide to the accuracy 
of the quadrupole polarizability. The first principles dipole polarizabilities are expected to 
be accurate to better than 0.5% and this has been confirmed by experiment. As discussed 
earlier, the uncertainty in the quadrupole polarizability for Na-like ions should be smaller 
than the dipole polarizability. Therefore it is not credible to postulate large errors in the 
atomic structure calculations of the quadrupole polarizability on the basis of a Bq derived 
from the polarization plot. It makes more sense to use the theoretical Cq to estimate the 
size of unaccounted systematic effects in the measured energy shifts. 

The large uncertainties in Bq do not detract greatly from the the accuracy of the dipole 
polarizability. One of the reasons higher-order effects can substantially impact Bq is that 
AEq is small because of the cancellation between a2 and (3i. However, the relatively small 
size of Bq means a large uncertainty in Bq has a relatively small impact on the derived ai. 

The impact of possible systematic errors on Mg"*" was determined by redoing the linear 
regression with a fixed value of Bq that was constrained to lie between —251.2 ± 19.0. The 
uncertainty of 19.0 was derived by adding the statistical uncertainty of 7.9 from the initial 
linear regression fit to |240.1 - 251. 2|, the difference between the Cq of Table |VTT1 and the 
initial Bq from the linear regression. This gave a revised uncertainty of 6B4 = 0.015, leading 
to a final ai of 35.04(3) a.u.. 
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The same analysis was repeated for Si^^. In this instance the derived value of ai was 
7.433(25) a.u.. 

F. Estimate of the resonant oscillator strengths. 

As the polarizabilities are dominated by the resonant transition it is possible to derive 



an estimate for the resonant multiplet strength 



801 1 ■ We use the relation 



C _ tt]^ ttcorc /„„x 

J3s-3p — 2 4 ■ V / 



In this expression ai is the polarizability extracted from the polarization plot while ttcore 
is the net core polarizability, and a'l is the valence polarizability excluding the resonant 
transition. For the Mg+ multiplet, we use 

35.044- 0.112- 0.463 

Sss-sp = : : = 8.439 . 30) 

3p 4.08436 ^ ^ 

Using the uncertainties detailed earlier, the final value is 8.439(11). This is equivalent to 
a line strength of S's^-sp^^^ = 11.25(2), in agreement with the recent experimental value of 
11.24(6) [50]. 

Repeating the analysis for Si^^ gave a multiplet strength of 3.519(16) for the 3s — > 3p 
transition. This is equivalent to S3s-3p.^/2 ~ 4.693(24) which is 0.14% larger than the MBPT- 
SD line strength of 4.686. 

VI. CONCLUSIONS 

A survey of polarization parameters of the Mg"*" and Si^^ ion states relevant to the 

nnn 

analysis of the RESIS experiments by the Lundeen group 15|, 16[ have been presented by 
two complementary approaches. The reanalysis of the fine-structure intervals gave dipole 
polarizabilities of 35.04(4) a.u. for Mg+ and 7.433(25) a.u. for Si^+ . The HFCP and MBPT- 
SD calculations give polarizabilities that lie within 0.2% of each other for Mg"*" and 0.3% 
for Si^+. The ab-initio MBPT-SD dipole polarizabilities of 35.05(12) and 7.419(16) a.u. 
respectively agree with the experimental dipole polarizabilities to accuracy of better than 
0.3%. 
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One notable feature of the present analysis is the very good agreement between the 
HFCP and MBPT-SD calculations. Indeed, the MBPT-SD calculation agrees better with 
the computationally simple HFCP calculation, than it does with two very large CI type 
calculations. For example, the polarizabilities of the completely ab-initio CI calculation 
[13I are about 1.5% larger than the MBPT-SD and HFCP polarizabilities. We conclude 
that a semi-empirical calculation based on a HF core can easily be superior to a pure CI 
calculation unless the CI calculation is of very large dimension. The HFCP approach has 
the advantage of tuning the model energy levels to experiment and this goes a long way to 
ensuring that many of the interesting observables will be predicted accurately. There is one 
feature common to the HFCP and MBPT-SD approaches. Both approaches approximate 
the physics of the dynamical corrections beyond HF/DF, but within those approximations 
an effectively exact calculation is made. 

There are two major sources of systematic error that can impact the interpretation of 
the RESIS experiment. To a certain extent one has to choose the {n,L) states to navigate 
between the Scylla [81] of non-adiabatic corrections and the Charybdis [81] of Stark shifts. 
If L is too small, then the AEj^s^sl shift becomes larger than AEq, thus invalidating the 
use of Eq. (Q. On the other hand. Stark shift corrections become increasingly bigger as L 
becomes larger. These problems are most severe in Si^"*" and are responsible for the slope 
of the polarization curve being different from the atomic structure predictions. An explicit 
two-state model of long range polarization interactions is probably needed to realize the full 
potential of the RESIS experiment 
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TABLE I: Theoretical and experimental energy levels (in Hartree) of some of the low-lying states 
of the Mg+ and Si^+ ions. The energies are given relative to the energy of the Mg^'^ and Sr^"*" 
cores. The experimental energies are taken from the National Institute of Standards and Technology 
database [43]. The HFCP energies should be interpreted as the J weighted average of the spin-orbit 
doublet. 



State 


Experiment 




MBPT-SD 


HFCP 






Mg+ 






3'Sl/2 


-0.552536 




-0.552522 


-0.552536 


3Pl/2 


-0.390015 




-0.390030 


-0.389737 


3P3/2 


-0.389597 




-0.389611 




4si/2 


-0.234481 




-0.234470 


-0.234323 


3c^5/2 


-0.226803 




-0.226772 


-0.226804 


3c?3/2 


-0.226799 




-0.226768 




4Pl/2 


-0.185206 




-0.185210 


-0.185014 


4^3/2 


-0.185067 




-0.185071 




4^5/2 


-0.127382 




-0.127374 


-0.127373 


4^3/2 


-0.127379 


Si3+ 


-0.127372 




3-S1/2 

J./ i 


-1.658930 




-1.658973 


-1.658928 


3pi/2 


-1.334120 




-1.334094 


-1.332738 


3p3/2 


-1.332019 




-1.331999 




/ z 


-0.928210 




-0.928138 


-0.928302 


3^3/2 


-0.928205 




-0.928134 




4si/2 


-0.775097 




-0.775104 


-0.774681 


4pi/2 


-0.664433 




-0.664421 


-0.663640 


4P3/2 


-0.663696 




-0.663684 




4^5/2 


-0.519810 




-0.519800 


-0.519743 


44/2 


-0.519809 




-0.519799 




4/5/2 


-0.501044 




-0.501044 


-0.501033 


4/7/2 


-0.501032 




-0.501035 
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TABLE II: Line strengths (in a.u.) for the resonance transitions of Na, Mg+, Al^+ and Si"^^. 
Experimental values with citations are also given. The MBPT-SD results for Na and Al^"*" are 
taken from flOl. 



Transition 



HFCP 



MBPT-SD 



BSR-CI 



Experiment 



S. 



(1) 

3s-3pi/2 

'3s-3p3/2 



S- 



(1) 

3s-3pi/2 

-"33-3^3/2 



s. 



(1) 

3s-3pi/2 



c(l) 
3s-3p3/2 



s. 



(1) 

3s-3pi/2 

:.(1) 
3s-3p3/2 



12.44 



24.88 



5.602 
11.20 



3.398 



6.796 



2.333 
4.666 



Na 



12.47 



24.94 



Mg-i 



5.612 
11.23 



A12+ 



3.404 



6.817 



Si 



^3+ 



2.341 
4.686 



12.60 



25.20 



5.644 
11.29 



3.422 



6.851 



2.350 
4.707 



12.412(16) [44, 45] 
12.435(41) [46] 
24.876(24) [47] 

24.818(34) [44, 45) 
24.844(54) [48] 



5.645(44) [49 
11.33(12) [49 
11.24(6) [50] 



3.01(29) [51] 
3.11(15) [52] 
3.31(35) [53J 
6.02(57) [51] 
6.35(45) [52] 



2.35(10) [54 
4.70(20) [54 
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TABLE III: Line strengths (in a.u.) for various transitions of Mg^ and Si'^'^. The Une strengths 
are mainly for dipole transitions with the exception of the 3s — > 3d and 3s Ad transitions. 







Mg+ 






Si^^+ 




Transition 


HFCP 


MBPT-SD 


BSR-CI 


HFCP 


MBPT-SD 


BSR-CI 


c-(l) 

•^3s-4pi/2 


0.00251 


0.00261 


0.00211 


0.0382 


0.0385 


0.0385 


0(1) 

'^3s-4p3/2 


0.00501 


0.00460 


0.00362 


0.0764 


0.0744 


0.0738 


■^38-5^1/2 


0.00395 


0.00402 


0.00366 


0.0138 


0.0139 


0.0138 


'^3s-5p3/2 


0.00790 


0.00763 


0.00692 


0.0276 


0.0270 


0.0268 


qW 

'^3pi/2-4s 


2.887 


2.868 


2.886 


0.6410 


0.6334 


0.6328 


0(1) 

•^3^3/2-48 


5.773 


5.779 


5.815 


1.282 


1.284 


1.283 


0(1) 
3pi/2-5s 




0.2117 


0.2115 




0.0633 


0.0629 


3p3/2-5s 




0.4247 


0.4243 




0.1284 


0.1267 


s 

3pi/2-3(i3/2 


17.32 


17.29 


17.35 


5.923 


5.933 


5.955 


0(1) 

3P3/2 -3^3/2 


3.463 


3.468 


3.482 


1.185 


1.190 


1.195 


3P3/2 -3^5/2 


31.17 


31.21 


31.33 


10.66 


10.71 


10.75 


^1 ^ 

S 

3Pl/2 -4(^3/2 


0.4100 


0.4168 


0.4631 


0.0234 


0.0212 


0.0204 


3p3/2 -4^3/2 


0.0820 


0.0825 


0.0918 


0.00468 


0.00455 


0.00438 


c(i) 

3^3/2-4^5/2 


0.7380 


0.7423 


0.8291 


0.0421 


0.0410 


0.0395 


^(2) 
3s— 3ci3/2 


97.52 


97.51 




16.81 


16.85 




0(2) 

'-^3s-3(i6/2 


146.3 


146.3 




25.21 


25.27 




q(2) 

'^3s-4(i3/2 


3.615 


3.638 




0.2732 


0.2659 




q(2) 

'^3s-4(iB/2 


5.422 


5.455 




0.4098 


0.3986 
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TABLE IV: The polarizabilities for the 3s, 3p and 3d states of Mg+ and Si^^. The tensor 
polarizabiUties are for the Mj = J states. For states with £ > 0, the MBPT-SD average values 
represent the weighted values for the spin-orbit doublet. 



State 


Qfl ( 


au) 




(au) 


Q!2 ( 


au) 




( au) 


HFCP MBPT-SD 


HFCP 


MBPT HFCP MBPT-SD 


HFCP 


MBPT-SD 




34.99 


35.05 


n 





156.1 


156.1 


1 71 5 


1 71 Q 

J- I J- 1/ 


iVlg \opi/2) 




31.60 


n 







340.2 




1 1 778 






31.88 


1.162 


1.156 




343.0 




11879 


Mg~'~(3p - Average) 


31.79 


31.79 






341.7 


342.1 


11839 


11845 


iVlg, V"-'"3/2 7 




189.3 


—78 47 


-79.15 




-9336 




2 857 X 1 0^ 






188.6 


-112.1 


—112.2 




—9341 




2.860x10^ 


Mg'*"(3d - Average) 


189.5 


188.9 






-9611 


-9339 


2.855x10^ 


2.859x10^ 


Si=^+(3s) 


7.399 


7.419 








12.13 


12.15 


47.03 


47.15 


Si'+(3pi/2) 




3.120 










13.05 




155.1 


Si-^+(3R3/2) 




3.183 


1.459 


1.462 




13.21 




157.1 


Si'^+ {3p - Average) 


3.158 


3.162 






13.17 


13.16 


156.3 


156.5 


Si3+(3d3/2) 




5.168 


-0.6083 


-0.631 




58.61 




695.2 


Si3+(3d5/2) 




5.131 


-0.8690 


-0.848 




58.61 




696.2 


Si^+(3(i - Average) 


5.135 


5.146 






58.43 


58.61 


693.2 


695.8 
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TABLE V: Breakdown of the different contributions to the dipole polarizabihties of Mg^ and Si 
The ep contribution includes both pseudo-state and continuum states. Dipole polarizabihties from 
other sources are also listed with citation. The estimated uncertainties for the different components 
of the uncertainty as estimated in brackets. The RESIS reanalysis are taken from the reanalysis of 
the RESIS fine-structure intervals described later. 



Quantity 




Mg+ 




Si3+ 


HFCP 


MBPT-SD 


HFCP 


MBPT-SD 


3s 3p 


34.413 


34.478(100) 


7.153 


7.180(6) 


3s U-6)p 


0.021 


0.020(0) 


0.054 


0.053(0) 


3s — > ep 


0.091 


0.087(4) 


0.030 


0.029(1) 


Core 


0.481 


0.481(10) 


0.162 


0.162(8) 


Core- Valence 


-0.018 


-0.018(2) 


-0.005 


-0.005(1) 


Total 


34.99 


35.05(12) 


7.394 


7.419(16) 


CI [13] 




35.66 




7.50 


RESIS [5, 6] 




35.00(5) 




7.426(12) 


Laser Exp [55] 




33.80(50) 






/-sums [56] 




35.1 






RCC [57] 




35.04 






RESIS reanalysis 




35.04(3) 




7.433(25) 



32 



TABLE VI: Breakdown of the different contributions to the quadrupole polarizabihties of Mg"*" 
and Si^+. The ep contribution includes both pseudo-state and continuum states. The quadrupole 
polarizabihties from a RESIS analysis is listed. 



Quantity 




Mg+ 




Si3+ 


HFCP 


MBPT-SD 


HFCP 


MBPT-SD 


3s 3d 


149.69 


149.68(32) 


11.502 


11.529(9) 


3s ^ (4-6)(i 


4.99 


5.01(4) 


0.240 


0.235(0) 


3s — > ed 


0.86 


0.85(6) 


0.289 


0.280(8) 


Core 


0.52 


0.52(6) 


0.102 


0.102(12) 


Total 


156.1 


156.1(5) 


12.13 


12.15(3) 


RESIS [5, 6] 




222(54) 






RCC [57] 




156.0 
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TABLE VII: The polarizabilities and C„ parameters computed from the composite list of HFCP 
and MBPT-SD matrix elements. The parameters tabulated here can be regarded as the recom- 
mended theoretical values. The C7, Cs and Csl parameters were used in the analysis the RESIS 

sp(-ctra for Mg^ and Si'^^. 



Quantity 


Mg+ 


Si3+ 


ai 


35.05(12) 


7.419(16) 


0:2 


156.1(5) 


12.15(3) 




1715(6) 


47.03(12) 


Pi 


106.0(3) 


11.04(1) 


P2 


236.1(5) 


8.065(6) 


71 


324.7(9) 


16.82(1) 


"112 


2416(52) 


89.74(41) 


"1111 


3511(90) 


51.19(28) 


C4 


17.53(6) 


3.710(8) 


Ce 


-240.1(12) 


-27.06(4) 


C7 


-1727(27) 


-125.6(2) 


Cs 


10672(92) 


553.1(6) 




1169(3) 


60.54(5) 



34 



TABLE VIII: Various energy corrections (in units of MHz) for the n = 17 intervals of Mg^ and 



the n = 29 intervals of Si'^"^ . These were computed using C„ values of Table IVIIl 



n 


Li 


L2 


AErel 




AEq 


AEj 


AEg 


AEsL 


AEscc 














Mg 


.+ 










17 


6 


7 


0.7314 


1555.7935 


-53.7751 - 


-20.4880 


7.1989 


31.5551 


8.1506 


-0.1122 


17 


7 


8 


5593 


678 8962 


—12 6733 


—3.4612 


0.8512 


4.9308 


1.5012 


—0.1702 


17 


8 


9 


0.4416 


326 8907 


—3.5557 


—0.7292 


0.1327 


0.9816 


3393 


—0 2320 


17 


9 


10 


0.3575 


169 8765 


— 1.1373 


—0.1808 


0253 


2325 


0896 


—0.2723 


17 


10 


11 


2953 


93 8138 


—0.4026 

Si^ 


— 0508 

+ 


0056 


0.0627 


0.0267 


—0 3039 


29 


8 


9 


7.2052 


1172.2322 


-67.1286 - 


-27.9320 


11.3800 


83.6138 


2.4220 


-0.1658 


29 


9 


10 


5.8328 


614.9317 


-22.2670 


-7.2876 


2.3124 


21.0913 


0.6574 


-0.3049 


29 


10 


11 


4.8184 


343.2123 


-8.2257 


-2.1728 


0.5526 


6.1258 


0.2026 


-0.4199 


29 


11 


12 


4.0474 


201.5328 


-3.3152 


-0.7211 


0.1502 


1.9896 


0.0693 


-0.6026 


29 


11 


13 


7.4953 


324.9889 


-4.7508 


-0.9824 


0.1956 


2.6972 


0.0951 


-1.3540 


29 


11 


14 


10.4675 


403.3750 


-5.4110 


-1.0843 


0.2106 


2.9685 


0.1054 


-2.3123 
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